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Abstract. We discuss recent results on the ionizing fluxes of O stars 
obtained from our "combined stellar structure and atmosphere models" 
{CoStar) accounting for stellar winds, non-LTE effects and line blanket- 
ing. The implications on the ionization structure of H ii regions are 
summarized, and observational constraints on the ionizing spectra and 
the total ionizing photon fluxes are presented. Using our CoStar models 
we derive new consistent predictions for the H and He i ionizing fluxes of 
steady-state massive star populations. Implications for the interpretation 
of observations of 'Diffuse Ionized Gas' in galaxies are discussed. Finally 
we present preliminary model calculations aiming to improve our current 
predictions further and compare our results to recent results from the 
Munich group. 



1. Introduction 

The "combined stellar structure and atmosphere models" (hereafter CoStar) 
models have been applied in detail to the main-sequence phase (Schaerer et 
al. 1996ab; paper I & II). Here we concentrate on some important predictions 
regarding the spectral evolution on the main sequence, namely recent results 
on the ionizing fluxes of OB stars for which considerable progress has been 
achieved in the recent years. Predictions of ionizing fluxes are of special interest 
for analysis of massive stars, H II regions, young starburst galaxies and related 
objects. 

We shall first present current models for O stars (Sect. ^). In Sect. ^ 
we will summarize some of the implications recent ionizing fluxes have on H II 
regions and how observations of H ii regions might be used to constrain model 
atmospheres. New predictions of ionizing fluxes of steady-state massive star 
populations and some of their implications on the diffuse ionized gas in galaxies 
are given in Sect. |4] In Sect. |5] we flnally present preliminary results from new 
calculations aiming to improve our current predictions further. 

2. Ionizing fluxes O stars 

The input physics of the CoStar models is described in detail in paper I. The 
most important ingredients regarding the prediction of ionizing fluxes are the 
treatment of non-LTE, the stellar wind, and line blanketing (see Schaerer & de 



Koter 1997, hereafter paper III). The current sets of models (paper III) cov- 
ers the entire main sequence evolution for stars with initial masses between 20 
and 120 Mq and includes a set for solar metallicity (Z=0.02) and low metal- 
licity (Z=0.004 and appropriately scaled wind parameters). The coverage ap- 
proximately corresponds to spectral types from 03 to BO and all luminosity 
classes. As discussed in paper III, reliable predictions for later types are cur- 
rently not available for physical and technical reasons. Spectral energy distri- 
butions, H II region models (cf. below), and additional information is a vailable 
on our Web page at http : / / www . stsci . edu/ ftp/ science/starburst and on 
CD-ROM (see Leitherer et al. 1996). 

Extensive comparisons of our predictions with previous calculations (Kurucz 
models, plane parallel line blanketed models) are given in paper III. A detailed 
discussion of the physical processes affecting the ionizing spectrum is also given 
there. The main results can be summarized as follows: 1 ) The flux in the He ii 
continuum is increased by 2 to 3 (3 to 6) orders of magnitudes compared to 
predictions from plane parallel non-LTE (LTE) model atmospheres (cf. Gabler 
et al. 1989). 2) The flux in the He i continuum is known to be increased due 
to non-LTE effects. In addition, the combined effect of a mass outflow and line 
blanketing leads to a flatter energy distribution in the He i continuum. This 
has profound implications on the structure of H ii regions (Sect. 3) To a 
lesser extent than at high energies, the flux in the Lyman continuum is also 
modified due to line blanketing and the presence of a stellar winds. Generally 
one obtains a reduction of the Lyman continuum flux compared to plane parallel 
model atmospheres. 

Based on our CoStar models we have also derived integrated photon fluxes 
for stellar parameters corresponding to the new spectral type and luminosity 
class calibration of Vacca et al. (1996). Our results are compared with former 
calibrations of Panagia (1973) and Vacca et al. (1996), both based on plane 
parallel model atmospheres. It must, however, be remembered that such a 
calibration is quite strongly dependent on the assignment between spectral types 
and stellar parameters, for which considerable uncertainties remain for O stars 
(cf. Vacca et al., Crowther 1997, Lamers et al. 1997, Schmutz, these proceedings). 



3. The impact of CoStar ionizing fluxes on H il regions and possible 
tests of atmosphere models 

Using the photoionization code PHOTO we have computed a large grid of single 
star H II region models based on the new CoStar ionizing fluxes (Stasihska & 
Schaerer 1997, hereafter paper IV). For comparisons a similar grid has been 
calculated using fluxes from plane parallel LTE model atmospheres of Kurucz 
(1991). The main additional parameter determining the structure of the H ii 
region, the ionization parameter U, is varied over a wide range corresponding 
to conditions expected between the most compact H ii regions and the most 
diluted interstellar medium. Detailed results are available on our Web page and 
on CD-ROM (see above). 

The flatter spectral energy distribution in the He i continuum (cf. Fig. |l| 
left panel) generally leads to higher ionic ratios in the nebula if CoStar models 
are compared to plane parallel atmospheres at the same stellar temperatures. In 
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Figure 1. Left panel: Comparison of the EUV continuous spectral 
energy distribution from a CoStar model (solid line) and a Kurucz 
model (dashed) for Tefj ~ 45000K (from paper IV). Right panel: Pre- 
dicted ionic Ne^"''/©^^ abundance ratio as a function of the 0"'"^/S^"'" 
ionic abundance ratio for the solar metallicity models (squares: CoStar 
models, triangles: Kurucz models). The solid lines connects photoion- 
ization models for a "standard" ionization parameter (cf. paper IV) 
and varying stellar temperatures. Dotted lines show variations of the 
ionization parameter. Observations from the IR line measurement of 
Simpson et al. (1995) are shown by filled circles. 



particular, with ionizing fluxes from CoStar models one obtains ionic fractions 
of x(N"^)/x(0"^) and x(Ne''"*')/x(0"'"^) which are close to unity for ionization 
parameters typical of H ii regions and over a large range of stellar temperatures. 
Therefore N^/O^ and Ne"^"*"/©"*"*" should provide safe indicators to derive N/0 
and Ne/0 abundance ratios over a wide range of astrophysical conditions. 

As shown by the exploratory calculations of Sellmaier et al. (1996) the use 
of ionizing fluxes from non-LTE wind models including line blocking seems to 
resolve the so-called [Ne iii] problem, which was made particularly clear when 
Simpson et al. (1995) plotted the values of Ne"''+/0+"'" as a function of 0"'""'"/S"'""'" 
(all derived from far-IR line measurements). The problem is illustrated in Fig. 
(right panel), where the observations indicate an essentially flat Ne"'"^/0"'""'" 
while the predictions of photoionizaton models using the Kurucz atmospheres 
fail to reproduce the observations. Indeed, the higher ionic ratios obtained from 
photoionization models using the CoStar fluxes are in much better agreement 
with the observations than models based on Kurucz atmospheres (see Fig. |^). 
As shown in paper IV a similar "observational test" can be performed by us- 
ing optical lines. Both comparisons support our ionizing fluxes from non-LTE 
atmospheres including stellar winds and line blanketing. 

Although such comparisons may represent extremely useful tests of pre- 
dicted EUV fluxes, good constraints on the ionization parameter and a precise 
knowledge about the ionizing sources (stellar type, number of stars) must be 
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obtained for accurate comparisons. This should be possible with future observa- 
tions of, in the ideal case, integrated single star H II regions. Careful studies of 
the metallicity dependence (and hence also the dependence on the wind proper- 
ties) of the predicted fluxes will be both necessary and extremely useful for many 
astrophysical applications. For example the great sensitivity of the Ne (and sev- 
eral other IR) lines on the shape of the ionizing spectrum in the He I continuum 
has important consequences for the interpretation of emission line spectra of H ii 
regions and young starbursts in the wavelength range observed by ISO. Derived 
stellar temperatures, IMF slopes and upper mass cut-offs depend strongly on 
the adopted ionizing fluxes as illustrated by Kunze et al. 1996. 

3.1. Constraints on the Lyman continuum production of OB stars 

Another fundamental test atmosphere models have to be subjected to is of course 
whether the integrated ionizing photon flux (e.g. in the Lyman continuum) is 
correctly predicted. The best example of such a test is certainly the B2II star 
e CMa, for which direct EUVE observations below the Lyman limit could be 
obtained (Cassinelli et al. 1995). For this case a large excess (by a factor of 30) 
of the Lyman continuum flux was found compared to plane parallel atmosphere 
models ! Although this excess can be explained by a temperature excess in 
the continuum forming layers (see Cassinelli et al., paper IV, Aufdenberg et al. 
1997) the physical cause for this temperature increase remains to be explained. 
Obviously in a larger context the fundamental question is, however, if the Lyman 
continuum excess is a general phenomenon, and if so for what type of stars ? 

In principle this question can easily be answered by observations measur- 
ing the integrated recombination rate (e.g. Ha luminosity) of nebulae around 
single stars of known spectral type. Detection of excess radiation compared to 
predictions for the determined spectral type could indicate an intrinsic excess 
of the Lyman continuum flux from the central star (or an undetected additional 
source). A lack of radiation can be interpreted as leakage of photons out of the 
region or important internal absorption by dust. To our knowledge no such sys- 
tematic comparison has been done so far, especially for B stars. Few statements 
can nevertheless already be made at this stage. 

For ultracompact H II regions which, according to Tcfj derived from FIR 
line ratios, contain stars of spectral types BO-07 a slight excess might be found 
for some cases (Afflerbach et al. 1997). However, their stellar content is not 
known directly and multiple sources can therefore not be excluded. The first 
ultracompact H II region where the spectral type of the central star has been 
directly identified (Watson k, Hanson 1997) indicates leakage or dust absorption. 
Oey (these proceedings) and Oey & Kennicutt (1997) have recently compared 
observed H II region luminosities with the Lyman continuum luminosities pre- 
dicted for their well determined stellar content. The vast majority is found to be 
leaking photons (up to ~ 50 %). From their samples, mostly dominated by sev- 
eral O stars, an excess of Lyman continuum flux can on average be excluded for 
O stars. Similar observations of less populated regions and nebulae dominated 
by later spectral types would be highly rewarding. 
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4. Photoionization of the diffuse ionized gas by massive stars 



The multi-phase interstellar medium (ISM) contains a diffuse ionized component 
representing a significant fraction of its total mass and volume. This gas is often 
referred to as the 'Reynolds Layer' in the Milky Way, or the 'Diffuse Ionized 
Medium' (DIM), 'Warm Ionized Medium' (WIM), 'Diffuse Ionized Gas' (DIG) 
etc. in external galaxies. We shall adopt the acronym 'DIG'. It is usually found 
that OB stars provide enough power to ionize the DIG (e.g. Reynolds 1984, 
Hunter & Gallagher 1990, Ferguson et al. 1996, Martin & Kennicutt 1997). 
The picture of photoionization by massive stars has, however, been challenged 
by recent measurements of a low He ionization fraction in the Galactic DIG 
(Heiles et al. 1996). On the other hand the observations of Martin & Kennicutt 
(1997) show that important variations in the physical state of the DIG exist 
among galaxies. 

To provide a quantitative base to address some of the questions raised by 
these studies it is highly desirable to use results from the most appropriate 
atmosphere and evolutionary models. For example, the expected He ionization 
of the DIG is dependent on the hardness of the ionizing fluxes of OB stars, which, 
as discussed above, has been revised considerably in recent years. Following the 
lines of Heiles et al. (1996) we will therefore derive new consistent calibrations of 
the H and He i ionizing fluxes for massive star populations based on our CoStar 
models. Some implications for the DIG of various galaxies will be discussed in 



Sect. 4.2 



4.1. Ionizing fluxes of massive star populations 

For studies of the DIG one usually considers a stellar population at a steady- 
state equilibrium ("constant star formation"). The basic quantities of interest 
are then (cf. Heiles et al. 1996): the lifetime r(m) of a star of initial mass m, its 
time-averaged Lyman continuum luminosity S}i{rn), and its time-averaged Hel 
continuum luminosity SHel("^)- These quantities can then be integrated over an 
initial mass function (IMF) to derive the total output of a stellar population. 

For the derivation of these quantities we use the CoStar models (papers 
I-III). We also extend the calculations to include the post main-sequence (MS) 
phases, which have been neglected in all previous calculations. The adopted 
input physics (tracks, atmospheres) is the same as in the synthesis models of 
Schaerer (1996) and Schaerer & Vacca (1997): we use the most recent Geneva 
tracks, Schmutz et al. (1992) atmospheres for Wolf-Rayet (WR) stars, and Ku- 
rucz (1991) atmospheres for the remaining phases. 

Fitting the results for solar metallicity models we obtain the following de- 
scriptions: 

logTMs("i) = 9.818 - 3.2861ogm + 0.8074(logm)^, (1) 

logrtotai("^) = 9.986 - 3.496 log m + 0.8942(log m)2. (2) 

tms and Ttotai (in years) are the main sequence and the total lifetime respectively. 
Eqs. H and ^ are valid for 7 < m/M© < 120. Note that the most massive stars 
(m > 85M0) enter the WR phase during core H-burning. In this case we have 
restricted the definition of the "MS phase" to the non-WR phase. We note that 



5 




Figure 2. Left panel: Hel/H hardness of the ionizing flux as a func- 
tion of initial mass. Time- averaged value over the entire lifetime (solid 
for solar metallicity; long-dashed for 1/5 Zq) and over the MS lifetime 
only (short-dashed for Zq) derived from our models. Dotted: Eq. 8 of 
Heiles et al. (1996). Right panel: Hel/H hardness of the ionizing flux of 
an integrated steady-state stellar population as a function of the upper 
mass cut-off. Solar metallicity models: Salpeter IMF (solid), a = 4. 
(dotted), a = 5. (short-dashed). Model at 1/5 Zq and Salpeter IMF: 
long-dashed. 



the MS lifetimes differ by less than ~ 15 % from the values of Glisten & Mezger 
(1982) used by Heiles et al. The mean Lyman continuum luminosity can be 
fitted by 

log su{m) = 32.47 + 17.35 log m - 4.362(log mf, (3) 

for the same mass range as above. It turns out that this expression represents 
well both the average over the MS and the total average. Our values agree well 
with those of Cox et al. (1986) for m ~ 20 — 40 Mq. Otherwise our values for sh 
are lower, with differences up to a factor of ~ 5 at the border of the fit interval. 

To describe the He i ionizing flux shcI we use C{m) = shci/sh (see Fig. ^), 
which can be approximated by 

((m) = -0.1604 + 0.0093m for 20 < m/M© < 40 (4) 
C{m)MS = 0.1656 + 0.0014m for 40 < m/M© < 120 (5) 
C(rn)totai = 0.1861 + 0.0007 m for 40 < m/M© < 120 (6) 

For m < 20 Mq the fraction of He i ionizing photons is negligible (C(W') < 0.01). 

The "time-averaged hardness" ^ as a function of the initial mass m is shown 
in Fig. ^ (left panel). Also shown is C(™') given by Heiles et al. (1996), which 
was obtained from the plane parallel non-LTE atmospheres of Kudritzki et al. 
(1991). Note that these atmosphere models predict a similar hardness (in terms 
of Qnei/Qu) of the ionizing flux as the CoStar models including also wind ef- 
fects and line blanketing. The data of Heiles et al., however, only correspond 



6 



to values typical for characteristic parameters of main sequence stars; any tem- 
poral evolution has been neglected. The evolutionary tracks in the considered 
mass range extend to temperatures below 40000 K on the MS. Given the rapid 
decrease of the He i ionizing luminosity below this temperature (e.g. paper III), 
the MS average of is systematically lower than the value given by Heiles et al. 
for m > 30Mq. For lower masses our results are in good agreement with the 
values of Heiles et al. Incidentally this is also the mass range which is critical 
for the interpretation of the observations of the DIG in the Milky Way (see Sect. 



1.2j ). It must, however, also be noted that this is likely the domain where cur- 
rent models are the most uncertain for the reasons discussed in paper III. Fig. |2| 
also illustrates the difference between averaging over the MS only and the more 
correct average over the MS and post-MS phases: For the most massive stars C 
averaged over their entire lifetime can be up to ~ 25 % lower than the MS value. 
This is due to their evolution in the WR phase, which leads to lower average 
luminosities. For comparison we have also plotted C (average over total lifetime) 
at 1/5 solar metallicity. As expected the average spectrum is slightly harder, 
mostly due to the temperature shift of the MS. In short we conclude that the 
time- averaged spectrum predicted from consistent modeling using in particular 
the ionizing fluxes from our CoStar models yields a lower Hel/H hardness than 
the results of Heiles et al. (1996)!^ 

Let us now consider the hardness C of the spectrum produced by an in- 
tegrated stellar population at a steady-state equilibrium. Figure § illustrates 
the dependence of C on the upper mass cut-off M-ap for various slopes a of the 
IMF. In our notation the Salpeter IMF is a = 2.35. Obviously the results are 
not sensitive to the lower mass cut-off Miow as long as Miow ~ 2OM0. Despite 
differences in all the ingredients (see above) our results agree well with those 
of Heiles et al. (their Table 6) for M^p < SOMq, whereas we obtain a softer 
spectrum for larger values of M^p. For illustration we also plot the run of C for a 
low metallicity population. We shall now briefly discuss some implications from 
these results. 



4.2. Implications for the diffuse ionized gas 

As shown by Heiles et al. (1996) the observations of the DIG in the Milky Way 
pose serious difficulties for the massive-star photoionization picture: while the 
high fraction of neutral He requires that very massive stars are excluded, the 
observed Galactic ionization requirement would imply a total star formation rate 
much larger than derived from other observations. Indeed the ratio of the Hel 
ionizing photon production rate Q(Hel) to (3(H) derived from the observations 
is < 0.013 - 0.027 (Heiles et al). 

Our results confirm entirely the interpretation of their Milky Way data. In 
contrast with some statements of Martin &; Kennicutt the use of our atmosphere 



^This contradicts several statements of Martin & Kennicutt (1997). Some comments must be 
made on the mass scales shown in their Fig. 3: 1) All scales are representative for dwarf 
stars and do not include the time average, whose importance was shown above. 2) Scale d was 
drawn from the preprint version of Heiles et al. which used LTE model atmospheres, while their 
final results are based on non-LTE models yielding a similar hardness as the CoStar models. 
Similarly for other results (Mup, Galactic star formation rate) the statements of Martin & 
Kennicutt refer to the preprint of Heiles et al., before important revisions were made. 
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models does not change the conclusions of Heiles et al. for the reasons discussed 
above (see footnote). Indeed, a comparison with the right panel of Fig. || shows 
that the observed Q(HeI)/(5(H) implies an upper mass cut-off of Afup ~ 25Mq. 
This result is robust even to large variations of the IMF slope (e.g. the very 
steep slopes found by Massey et al. (1995) for field stars) as also show by Fig. |2|. 
For a = 2.35 we also obtain the same Lyman continuum production as Heiles 
et al., which confirms the requirement of a large Galactic star formation rate to 
reproduce the observed Galactic ionizing photon production rate. 

As an unconventional explanation for these problems Heiles et al. consider 
the possible excess of Lyman continuum flux from B stars, motivated by the 
observed excess in the B2II star e CMa. If common among B stars this could 
potentially solve the above mentioned difficulties. As discussed in Sect. 3.1.| it 
should be possible to constrain the ionizing fluxes of B stars observationally in 
the future. 

The observations of DIG of three Magellanic irregular galaxies of Martin & 
Kennicutt (1997) reveal a different picture. There the He ionization fraction is 
approximately equal to the hydrogen fraction, which translates to a lower limit 
of Q(Hel)/Q(R) > 0.1 (Martin & Kennicutt). From the right panel in Fig. ^ this 
translates to a lower limit for the upper mass cut-off of Mup ^ 30-40 Mq for a 
Salpeter IMF. We also note that the high He ionization in these objects could in 
principle even be explained by a population with an IMF slope as steep as those 
derived by Massey et al. (1995) for field stars. For this case the required upper 
mass cut-off is Mup > 40-60 Mq, well compatible with numerous independent 
observational constraints (cf. Leitherer 1997, Massey 1997). The ionization by 
radiation leaking out from H ii regions (with a "normal" IMF slope) is, however, 
a more likely explanation, as shown by recent observations (Ferguson et al. 1996, 
Martin & Kennicutt 1997, Oey & Kennicutt 1997, Wang et al. 1997). 



5. Reliable ionizing fluxes for OB stars: further developments 

To improve our atmosphere models further and to compare the predicted spectra 
with UV observations we have recently upgraded the non-LTE calculations by 
including a detailed treatment of CNO and Si in addition to H and He (following 
de Koter et al. 1997). Improvements are also brought about for the completeness 
of the line list used in the line blanketing Monte-Carlo simulations. In this 
Section we shall briefly present some preliminary results on the EUV fluxes 
and compare our latest calculations with new models from the Munich group 
(Pauldrach, these proceedings). 

The main changes with respect to paper IV are the explicit non-LTE treat- 
ment of CNO and Si, whose populations were previously described by our mod- 
ified nebular approximation (see paper I and references therein). By the same 
token this also allows us to derive a finer description of the radiation temperature 
as a function of the ionization potential, which is used to describe the ionization 
structure of the remaining elements. The bound-free opacities of CNO and Si 
are now also included. In Fig. ^ we plot a comparison for two models from 
paper IV. For A ^ 310 A the new calculations agree well with the previous ones. 
Between 310 A and the He ii edge (228 A) our new models predict less flux, 
mostly due to the inclusion of the bound-free CNO opacity. Predictions for the 
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Figure 3. Comparison of the emergent fluxes from different models. 
Long-dashed: from paper IV, Sohd: new model (continuous distri- 
bution: thick, detailed line-|-cont.: thin). Dotted: Pauldrach model 
(smoothed spectrum traced by thin solid line), Dashed-dotted: Kurucz 
model. Left panel: model C2 (Tgfj ~ 40000 K), right panel: model El 
(Teff ~ 50000 K), 



He II continuum remain subject to other considerable uncertainties (see paper 
IV). 

In Fig. ^ we have also plotted the flux from two new atmosphere models 
including line blocking calculated with similar parameters (see Pauldrach, these 
proceedings). The overall agreement between the models is quite good for A > 
310 A, although some distinct features exist. At shorter wavelengths, however, 
the flux from the Munich model decreases more strongly, remaining fairly close to 
the Kurucz model for T^g ~ 40000 K. (left panel). We have done test calculations 
of photoionization models using both sets of fluxes. Regarding the nebular 
Ne"*"^/©^^ abundance ratio (cf. Sect. |3]) we obtained only minor changes for 
our models, while the Pauldrach model at 40000 K yield Ne++/0++ 0.1-0.15 
and 0^~'"/S~'"^ ^ 25, more similar to Kurucz models. 

Given important differences between the model techniques and the some of 
the input physics (e.g. line lists) it is difficult to explain the reason(s) for the 
differences in the predicted spectra. More detailed work remains to be done to 
further explore the dependence on various model assumptions and to improve 
the treatment of line blanketing and non-LTE effects in spherically expanding 
atmospheres. Concomitantly it appears very useful to pursue further studies of 
"simple" H ii regions to constrain atmosphere models as much as possible. In 
view of the difficulties hot star models have faced so far in reproducing correctly 
the observable UV from consistent calculations (see e.g. Groenewegen & Lamers 
1991, Pauldrach et al. 1994) such new tests could be crucial to reach the goal of 
obtaining accurate predictions of the spectral evolution of early type stars. 
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